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ABSTRACT 



We made VLBI observations at 8.4 GHz between 1997 and 2005 to estimate 
the coordinates of the "core" component of the superluminal quasar, 3C 454.3, 
^ ■ the ultimate reference point in the distant universe for the NASA/Stanford Gyro- 

scope Relativity Mission, Gravity Probe B. These coordinates are determined rel- 
ative to those of the brightness peaks of two other compact extragalactic sources, 
'ES ■ B2250-I-194 and B2252-I-172, nearby on the sky, and within a celestial reference 

frame (CRF), defined by a large suite of compact extragalactic radio sources, and 
nearly identical to the International Celestial Reference Frame 2 (ICRF2). We 
^ ■ find that B2250-I-194 and B2252-I-172 are stationary relative to each other, and 

' also in the CRF, to within la upper limits of 15 and 30 yuas yr~^ in a and 5, 

^ ■ respectively. The core of 3C 454.3 appears to jitter in its position along the jet 

■ direction over ~ 0.2 mas, likely due to activity close to the putative supermas- 

O ' sive black hole nearby, but on average is stationary in the CRF within la upper 

limits on its proper motion of 39 /las yr~^ (1.0c) and 30 /las yr~^ (0.8c) in a and 
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5, respectively, for the period 2002 — 2005. Our corresponding limit over the 
longer interval, 1998 — 2005, of more importance to GP-B, is 46 and 56 /xas yr~^ 
in a and 5, respectively. Some of 3C 454. 3's jet components show significantly 
superluminal motion with speeds of up to ~200 /las yr~^ or 5c in the CRF. The 
core of 3C 454.3 thus provides for Gravity Probe B a sufficiently stable reference 
in the distant universe. 

Subject headings: binaries: close — radio continuum: stars — stars: activity — 
stars: imaging — stars: individual (IM Pegasi) — techniques: interferometric 

1. Introduction 

Gravity Probe B {GP-B) is the spaceborne relativity experiment developed by NASA 
and Stanford University to test two predictions of general relativity. The experiment used 
four superconducting gyroscopes, contained in a low-altitude, polar orbiting spacecraft, to 
measure the geodetic effect and the much smaller frame-dragging effect. According to general 
relativity, each of these effects induces precessions of the gyroscopes in planes perpendicular 
to each other. For the geodetic effect, which depends directly on the Earth's mass, the 
predicted precession is 6.6 arcsec yr^^ and for the frame-dragging effect, which depends 
directly on the angular momentum of the Earth, it is 39 mas yr~^. GP-B was expected to 
measure each precession with a standard error <0.5 mas yr~^ relative to the distant universe. 
Because of technical hmitations, the spacecraft could not measure the precessions directly 
relative to the distant universe but only to an optically bright star, the guide star, chosen 
to be IM Pegasi (HR 8703). We must therefore determine IM Peg's proper motion relative 
to the distant universe, which is, for our purposes, best represented by extragalactic radio 
sources. 

For our part of the GP-B project, we determined the coordinates and the proper motion 
of the guide star in the radio relative to the "core" of the quasar 3C 454.3 (B2251-I-158). 
This core was tied to two other radio sources, which are compact, extragalactic, nearby to 
it on the sky, and also tied to a celestial reference frame (CRF) defined by a large suite 
of extragalactic sources. These ties are the main subject of this paper. Most important 
for GP-B, of course, is the bound that we place on the proper motion of the core, which 
serves as the principal reference for determining the proper motion of IM Peg. The GP-B 
project needs the proper motion of the optical source in IM Peg. The radio source in IM Peg, 
however, moves erratically with respect to the optical source. In order to be able to average 
as well as feasible over the erratic motion, we place our VLBI limit on the motion of the 
core of 3C 454.3 over as long a period as feasible. By contrast, for astrophysical purposes. 
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we place a more stringent bound on the core's motion, but only for a substantially shorter 
period of time. 

Apart from its relevance for GP-B, our observations and astrometric analysis are also of 



astr ophysical interest. 



Pauliny-Toth et al. 



he qu asar 3C 454.3 is a highly active superluminal radio source (see, 



198 7). It c onsists of a relatively compact region from which a bent jet 



e.g.: 

emanates (e.g. . lPauliny-Tothlll998[ l . Superluminal motion refers to those apparent transverse 
velocities of the components within the source that are measured to be greater than c, the 
speed of light. On the basis of synchrotron radiation theory, the core is generally identified 
as that component which is compact and has a flat or inverted radio spectrum. Observations 
at 43 and 86 GHz show tha t 3C 454.3 has suc h a component located in the eastern part of 
the brightness distribution ( iPagels et al.l 120041 ). Its characteristics are consistent with those 
expected for the environment of a supermassive black hole or the base of an associated jet, 
but not conclusively diagnostic of either one. An additionally powerful probe for the location 
of a possible supermassive black hole in an extragalactic source is to identify the component 
in the radio structure of the source which shows the smallest motion of all components or 
is stationary on the sky. Such a component would be a strong candidate for being closely 
related to the purported supermassive black hole which is likely to be both close to the center 
of mass of the source and virtually stationary on the sky. 

Placement of stringent limits on the proper motions of quasars and other compact extra- 
galactic radio sources distributed across the sky are being made by others on a routine basis 
through astrometric/geodetic VLBI observation^. Random errors of frequently observed 
sources may be as low as 6 /las, although systematic errors, mainly due to unaccounted 
propagation effects and source structure are believed to be in the range of 50-1000 /zas. 

All geodetic VLBI measurements are based on interferometric group delays. More ac- 
curate measurements can be made, for example, with targeted differ ential VLBI observa - 
tions of two or more sources located close to each other on the sky (IShapiro et al.lll979l ). 
In such cases, interferometric phases or phase delays can be used, yielding relative posi- 
tions and proper motions of properly selected celestial sources wit 
10 jUa s yr~^, respectively (e.g . , 
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Rioia &: Porcad l2000l : iBietenholz et al.l l200ll : iFomalont fc Kopeikin 
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Bartel et al. 



2003 
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1986 



Brunthaler et al. 



This paper is the third in a series of seven papers reporting on the astrometric support 
for GP-B for the purpose of defining the cosmological reference frame for the gyroscope 



^Apart from source positions, these observations yield antenna coordinates and velocities, and series of 
Earth orientation parameters. 
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precession measurements. I n the first paper of t his series we gave an overview of the astro- 
nomical support for GP-B (jShapiro et al.ll201ll . Paper I). In the second paper we focused 
on the characteristics of quasar 3C 454.3 and the other two extragalactic reference sources, 
B2250-I-194 and B22 52+172, and reporte d on their structure and structure changes with 
time and frequency (IRansom et al.ll2011al . Paper II). In this paper (Paper III), we report 
on the degree of stationarity of the core of the quasar 3C 454.3, which is the reference 
source for the guide star IM Peg and therefore pivotal for GP-B. In the fourth paper we 
present our VLBI a strometry analysis technique and compare it with other such techniques 
( iLebach et al.ll201ll . Paper IV). In the fifth paper we presei it our results for th e proper mo- 
tion and parallax of IM Peg relative to the core of 3C 454.3 (IRatner et al.ll201ll . Paper V). In 
the sixth p aper we report on the orbital motion of IM Peg and interpret the radio structure 
of the star ( IRansom et al.ll2011bl . Paper VI). Finally, in the seventh paper, we focus on the 
indiv idual epochs of observ ation of IM Peg and include a movie of the radio images of this 
star ( iBietenholz et al.ll201ll . Paper VII). 



Here we first briefly describe our observations, in § [2l We give characteristics and show 
representative images of 3C 454.3, B2250+194, and B2252-M72 in § 13 We describe our 
astrometry program in § |H We present astrometric results in § El El and [71 discuss these 
results in § [HI and give our conclusions in § [9] 



2. Observations 

As one of the strongest-emitting quasars at radio frequencies, 3C 454.3 has been observed 
in geodetic group-delay VLBI sessions since 1979. For our GP-B VLBI program we made 
use of observations from the total of 1119 such sessions between 1980 and 2008. In addition 
we used geodetic observations of B2250+194 from a total of 38 sessions between 1996 and 
2008 that were made in support of GP-B VLBI. 

The bulk of our GP-B VLBI efforts were devoted to phase-delay VLBI observations of 
IM Peg and our three reference sources, 3C 454.3, B2250+194, and B2252-I-172. A detailed 
description of these latter observations was given in Paper II; however, for the convenience 
of the reader, we give a summary here. 

We obtained 35 sets of 8.4 GHz VLBI observations in support of GP-B between 1997 
January 16 and 2005 July 16. We used a global array of 12 to 16 radio telescopes, which most 
often included MPIfR's 100 m telescope at Effelsberg, Germany; NASA/Caltech/JPL's 70m 
DSN telescopes at Robledo, Spain, Goldstone, CA, and Tidbinbilla, Australia; NRAO's ten 
25 m telescopes of the VLBA, across the U.S. A; NRAO's phased VLA, equivalent to a 130 m 
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telescope, near Socorro, NM; and, at early times, NRCan's 46 m Algonquin Radio Telescope 
near Pembroke, ONT, Canada, and, at later times, NRAO's 110 m GBT in WV. In each 
session we made interleaved observations of 3C 454.3, IM Peg, and B2250+194 by using a 
sequence of typically 3C 454.3 (80 s) - IM Peg (170 s) - 2250 (80 s). For the last 12 sessions, 
starting 2002 November, we also observed B2252+172, but only after every second sequence 
to allow greater concentration on the main three sources. The new observing sequence was 
3C 454.3 (80 s) - IM Peg (125 s) - 2250 (80 s) - 3C 454.3 (80 s) - IM Peg (125 s) - 2250 (80 s) 
- 2252 (90 sj^. In three sessions we also observed at 5.0 and 15.0 GHz. All observations were 
recorded in both right and left circular polarizations and processed on the VLBA hardware 
correlator at Socorro. 



3. The Celestial Reference Sources 

3.1. Sky Positions and Cosmological Distances 

In Figured] we show the positions of 3C 454.3, B2250+194, and B2252+172 along with 
that of the GP-B guide star, IM Peg. All four sources are located approximately along 
a single north-south axis allowing us to make easier use of interpolation to estimate and 
potentially reduce the contributions of the troposphere and the ionosphere to the total error 
in determining the sources' relative positions. In Table [1] we give the sky separations of 
the two extragalactic sources from 3C 454.3, the sources' flux densities, redshifts (when 
known), and angular diameter distances, the latter assuming an inhomogeneous Friedmann- 
Lemaitre-Robertson- Walker cosmology. For comparison, we also give the characteristics of 
IM Peg. 



3.2. Representative Images 



In Figures HI El and S] we show representative images of 3C 454.3, B2250+194, and 
B2252+172. The source 3C 454.3 is a superluminal quasar with the highest radio flux 
density of the three sources. It consists of a core region, which is primarily extended east- 
west, and can be well modeled for each of our 35 epochs by two compact components, CI and 
C2, separated by ~0.6 mas. When stud ied at a higher resolution, e.g., at 43 and 86 GHz, 



the same region has a complex structure (IGomez et al.lll999t iJorstad et aL 



12001 



2005), with 



CI being essentially unresolved and having a size at 86 GHz of <70 /ias ( jPagels et al.ll2004j ). 



^Here and hereafter we sometimes use as abbreviations 2250 for B2250+194 and 2252 for B2252+172. 
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Fig. 1. — A sky chart with coordinates of the three reference sources and the guide star, 
IM peg. The hnear scale is the same for right ascension (a) and declination (5) for the 
declination of 3C 454.3. 
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C2 is extended. Also, between these two component s there are others that ap pear to move 
away from CI toward C2 with superluminal speeds ( iJorstad et al.ll200ll . |2005| ): these other 
components are not individually visible in our lower-resolution 8.4 GHz images. 

Further to the west there are components, Dl and D2, that are visible only at later 
epochs in our observations. We modeled these as compact components. They bridge the 
core region to the outer part of the jet with a more than 10 mas extent in our images. 
This outer part of the jet bends toward the northwest and in fact ex t ends as far as a few 
arcseconds away from the core (see also ICawthorne fc Gabuzdal Il996l : iPauliny-TothI Il998l ) . 
The brightness peak of the extended 10 mas jet in our image is clearly visible at each of our 
35 epochs, and we modeled it as a compact component which we call Jl. The extended jet 
can be modeled as an elliptical Gaussian, which we call Jext. 

The source B2250+194 is 20 times weaker in terms of flux density than 3C 454.3, but 
much more compact, consisting of a central component with north-south extensions and an 
apparent length of ~5 mas. The source B2252+172 is the weakest of the three in terms of 
flux density but also the most compact, consisting of one dominant component and a barely 
visible extension to the west. For more detail on the structure of these sources and their 
evolution or lack thereof, see Paper II. 



3.3. Selection of 3C 454.3 as a Reference Source for GP-B 



Because of these differences in characteristics and separations from IM Peg, each of 
these sources has advantages and disadvantag reference source for IM Peg. The 

quasar 3C 454.3 has the advantage of the highest flux density and closest proximity to the 
guide star. The latter point is of chief importance, since the standard errors of astromet- 
ric VLBI measurements scale approximately linearly with the separa tion of the reference 
source from the target source (I Shapiro et al.lll979l : iPradel et al.ll2006l ). Another advantage 
i s that 3C 454.3 was used as the astrometric reference source for IM Peg as early as 1991 
( iLestrade et al.lll999l ). thus extending the time baseline of VLBI position determinations and 
allowing increased accuracy in our proper-motion estimates. The main disadvantage of this 
reference source is its complex structure. 

The advantage of the other two sources is their compact structure. However, B2250+194 
is relatively far away on the sky from IM Peg, and B2252-I-172 has a relatively low flux density. 
All things considered, we decided to use 3C 454.3 as the primary reference source for IM 
Peg. Our task was then to find the component in 3C 454.3 that is most closely associated 
with the dynamical center of the quasar and to test the stationarity of this component with 
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Aa (mas) 



Fig. 2. — An image of 3C 454.3 from observations on 2005 May 28 with components CI, C2, 
Dl, D2, Jl, and Jext indicated. The contours start at 10 mJy bm^^ and increase by factors 
of towards the peak. The peak brightness is 2.68 Jy bm~^ The rms brightness of the 
background noise is 0.73 mJy bm~^. The full- width at half- maximum (FWHM) contour of 
the Gaussian convolving beam is given in the lower right. North is up and east to the left. 
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Fig. 3. — An image of B2250+194 from observations on 2005 January 15. The contours are 
at 0.3, 0.5, 1, 2, 5, 10, 20, 40, and 80% of the peak brightness of 0.43 Jy bm'^ The rms 
brightness of the background noise in 0.08 mJy bm~^. The FWHM contour of the Gaussian 
convolving beam is given in the lower left. North is up and east to the left. 
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Fig. 4. — An image of B2252+172 from observations on 2005 January 15. The contours are 
at 3, 5, 10, 20, 30, ... 90% of the peak brightness of 0.012 Jy bm~^. The rms brightness of 
the background noise is 0.05 mJy bm~^. The FWHM contour of the Gaussian convolving 
beam is given in the lower left. North is up and east to the left. 
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respect to our two other reference sources and our CRF. 



Images of 3C 454.3 at 8.4 and 15 GHz (Paper II) and at 43 and 86 GHz flPagels et al. 



20041 ) show that the easternmost component, CI, remains compact at the highest frequencies 



and angular resolutions yet investigated, and has a flat or inverted spectrum in this frequency 
region. Other components or condensations show structure at 86 GHz, and, in cases where 
it was determinable, a steep spectral index. These characteristics indicate that for our 8.4 
GHz images, CI is likely the component most closely related to the putative supermassive 
black hole and the quasar's center of mass. 



4. The Astrometry Program for GP-B: Goal, Strategy, and Procedure 

Goal: Our main goal is to determine Cl's position and especially a bound on its proper 
motion relative to the distant universe, to confirm its suitability as the primary reference for 
GP-B. 

Strategy: The distant universe is for our purposes most usefully represented by compact, 
extragalactic radio sources. We choose a procedure where we first determine the position 
and proper motion of CI relative to our two reference sources, B2250+194 and B2252+172, 
and second relative to our CRF. For the first step, interferometric phase-delays are used 
exclusively. This has the advantage of simplicity, utilization of the same type of data (phase- 
delays) and analysis technique, and highest precision and possibly highest accuracy of the 
results. For the second, the results from the first are added to position and proper-motion 
determinations of B2250+194 in the CRF based on interferometric group delays from the 
geodetic VLBI sessions. This has the advantage of having as a reference not only two very 
distant and compact sources, but ~4000 sources that define our CRF. 

Procedure: We determine the position and the limit on proper motion of: 

1. CI relative to B2250+194 and B2252+172, and B2252+172 relative to B2250+194 
in the senses (CI - 2250), (CI - 2252), and (2252 - 2250) using our interleaved 
phase-delay VLBI observations. The combined result for CI is obtained as a weighted 
mean of the first two differenced solutions, with the third serving to demonstrate the 
consistency of our results and zero proper motion within the errors for the two sources 
relative to each other. Due to the source's compactness, any motions or brightness 
distribution changes of these two reference sources would likely be very small and 
therefore have only a marginal, if any, effect on our astrometric results for CI. 
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2. B2250+194 in the CRF using routine geodetic and astrometric group-delay VLBI ob- 
servations distributed by the International VLBI Service for Geodesy and Astrometr}!^!; 

3. B2252+172 in the CRF by adding the result from 2. to that from 1. in the sense (2252 
— 2250) + 2250, and confirming the position result and proper motion limit by using 
recent single-epoch geodetic group-delay VLBI observations of B2252+172; and 

4. CI in the CRF in two ways: First by adding the results from 1. to those from 2.: CI 
= (CI — 2250) + 2250, and second by adding the results from 1. to those of 3.: CI = 
(CI — 2252) + 2252. The combined result is obtained as a weighted mean from these 
two ways. 



5. Astrometric Results (1): Position Determinations for Each Observing 

Session 



5.1. Analysis of Interleaved Phase-delay Observations of 3C 454.3, 

B2250+194, and B2252+172 

Our VLBI data for the reference sources, 3C 454.3, B2250-M94, B2252+172, and also 
from IM Peg, were analyzed with an astrometric software package that was developed specif- 
ically for the analysis of the GP-B VLBI data. It includes a phase-connection program that 
automatically resolves 2tt ambiguities that exist in the set of VLBI phases for each baseline 
so as to c onvert them to pha s e delays. The sof tware also includes in the phase-delay fit- 
ting (e.g., IShapiro et al.lll979l : iBartel et al.l 119861 ) a Kalman filter (see Paper IV) to model 
the variations of th e troposphere, the ionosphere, and the clock offsets at each VLBI site 
(iLebach et al.lll999l ). In addition, we used two different models to initially correct for the 
effects of the ionosphere, one we call "JPL" which is part of NRAO's imaging package, AIPS, 
and based on GPS data provided by JPL, and the other, the older pa rametrized ionos pheric 
model "PIM," developed at USAF Research Labs (described by, e.g.. ICampbelllll999l ). 



To relate the phase delays to a particular reference point in 3C 454.3, namely the 
core component CI, all phase delays from differential astrometry involving 3C 454.3 were 
corrected for the structure of 3C 454.3, as represented by the CLEAN components produced 
with AIPS; CI served as the phase reference point. The other two extragalactic sources were 
deemed sufficiently compact for our purposes so that the reference point for each could be 



■^Available at 



http://ivscc.gsfc.nasa.gov/products-data/index.htinl 
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assumed to be the brightness peak in its image. We elaborate on this method of astrometric 
VLBI data analysis and compare it to other methods in Paper IV. 

5.2. Positions of the Components of 3C 454.3 Relative to Those of 

B2250+194 and B2252+172 

To test the positional stability of 3C 454. 3's component CI relative to our two reference 
sources, we determined, for each of our 35 sessions of 8.4 GHz observations, the coordinates 
of CI, and for comparison also those of C2, Dl, D2, Jl, and Jext, all relative to the brightness 
peak of B2250+194 and, for the last 12 epochs, also to that of B2252+172. We obtained two 
sets of coordinates by correcting for the effects of the ionosphere in two different ways, one 
set by using the JPL model and another set by using PIM. Although the JPL model was 
not available for our first 8 epochs but only from 1998 September 17 onward, we nevertheless 
in this paper use mainly phase-delay data corrected with the JPL model since it proved to 
be superior in that it resulted in smaller uncertainties of our astrometric estimates despite 
precluding the use of our earliest phase-delay data. We elaborate on the comparison below. 
(In Paper V we use PIM instead, because in that paper errors in modeling the ionosphere 
play a less significant role than they play here, and because PIM has the advantage that it 
can be used for all of our VLBI data.) We list our coordinate determinations for CI, C2, 
Dl, D2, Jl, and Jext relative to the brightness peak of B2250+194 in Table [2] and to that 
of B2252+172 in Table |3l all obtained with the JPL model. 

5.3. The Uncertainties of the Relative Positions of the Components of 3C 

454.3 

The uncertainties of the coordinates listed in Tables |5] and [3] were determined partly 
empirically, namely by adding a constant in quadrature to the statistical standard errors so 
as to obtain a reduced Chi-square of unity ix^/^ = where u is the number of degrees of 
freedom) in our residuals after solving for relative position and proper motion in a and 6 
separately. This constant is assumed to approximately reflect non-statistical errors. Accu- 
rate standard errors are difficult to estimate in any other way. They contain contributions 
from noise and from systematic errors, with the latter due mostly to deconvolution, source 
structure, and atmospheric and ionospheric variations. We next discuss and approximately 
quantify each contribution in turn. 
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5.3.1. Noise 

Noise in an image has a number of sources. The rms background brightness in the images 
is dominated by contributions from statistical noise in the radio signals and thermal noise in 
the receivers. However, for our relatively high dynamic range images (typically over 2,000 to 
1; see, e.g.. Figures |2]to|l]), the various uncertainties in the bright parts of the image are larger 
than their corresponding rms background brightnesses and are dominated by contributions 
which are not strictly random such as residual calibration errors and deconvolution errors. 
Given the small rms of the background brightness relative to the peaks in the maps, we 
conclude that this noise causes errors in the estimate of the separation of each of components 
CI and C2 in 3C 454.3 from B2250+194 by a correspondingly small portion of the HWHM 
(half-width at half-maximum) of the beam, namely by <5 //as, and of each of components 
Dl, D2, Jl and Jext in 3C 454.3 from B2250+194 by <10 /xas. The corresponding errors 
in our estimates of the separation of these components from the weaker source, B2252-(-172, 
are dominated by the source's lower peak-to- noise ratio but are still <10 /xas. 



5.3.2. Deconvolution errors 

Deconvolution errors are caused by the visibility measurements not filling the u-v plane 
of the VLBI array up to its highest angular resolution and by the resulting generation of side 
lobes in the image plane, which are not completely eliminated through the deconvolution 
process. We studied this type of error by using a noise-free model image similar to the image 
of 3C 454.3 at 8.4 GHz, Fourier-transforming the model to the u-v plane, and then using 
the same w-w-plane sampling as in one of our typical observing sessions. The generated 
u-v model data were then used for imaging and deconvolution. We then determined the 
difference between the position of each component in the model image and the position of 
the corresponding component in the deconvolved image. We found that the deconvolution 
error for each coordinate of each of the six components of 3C 454.3 is typically 30 /las 
and never larger than ~ 40 ixas. The deconvolution errors for sources like B2250+194 
and B2252+172 with relatively simple brightness distributions are doubtless smaller given 
the same u-v coverage as for 3C 454.3. We therefore conclude that the standard error in 
each coordinate of the separation of any component in 3C 454.3 from either B2250-I-194 or 
B2252+172 is typically 30 //as. 
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5.3.3. Structure errors 



Structure errors are caused by a mismatch between the Gaussian component model and 
the brightness distribution of the source. If the source were completely unresolved, a fit of a 
Gaussian to the image with the parameters of the convolving beam would give the position 
of the source with an essentially zero structure error. Indeed, the sources B2250+194 and 
B2252+172 are rather compact, with brightness peaks that can be clearly identified and 
located, and hence we expect the structure errors to be small in comparison to those for the 
components in 3C 454.3. As can be seen in Figure |2l components CI and C2 are located close 
together in comparison to the size of the beam and can barely be distinguished at several 
epochs (see Paper II). In addition, there is "confusing" emission in the neighborhood of 
these two components. These two characteristics cause significant structure errors, difficult 
to estimate quantitatively, but likely as large as a good fraction of the HWHM of the beam. 
The components Dl and D2 are weaker than CI and C2, and visible only at later epochs, 
but then clearly distinguishable. Their identification is likely as uncertain as that of CI 
and C2. Component Jl is always visible as a single peak; however, it is embedded in the 
extended brightness distribution of the 10 mas long jet rendering the identification of a 
component in that jet as uncertain as, for many epochs, the identification of CI and C2. 
The identification of Jext is more uncertain than that of any of the other components since 
this component represents the extended brightness distribution that Jl is embedded in. The 
extension is largely toward the northwest, thus affecting both position coordinates of Jext. 
A detailed analysis of these uncertainties is difficult to carry out quantitatively, in view of 
the unknown characteristics of the relevant structures. Based on our experience, however, 
we expect structure errors in each of the position coordinates to be ~ 30 /ias for components 
CI, C2, Dl, D2, and Jl and ~ 50 /xas for the component Jext. In comparison, the structure 
errors for B2250+194 and B2252+172 are negligible. 



5.3.4- Residual propagation-medium errors 



Using our Kalman filter (ILebach et al.ll201ll Paper IV) removes a large portion of the 
propagation-medium errors from the estimates of separation between two sources. Neverthe- 
less, errors in modeling the troposphere and the ionosphere at each site likely still represent 
the largest sources of error in our estimates of relat ive positions. These errors scale approx - 



imately with the separation between the sources (IShapiro et al.l Il979l : iPradel et al. 



20061). 



These errors also depend on the model used for the ionosphere (JPL or PIM). 



In our case the sources are oriented approximately north-south with the separations in 
a of B2250+194 and B2252+172 from CI being about equal, but with the separation in 6 
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of B2252+172 from CI being only about 40% of that for B2250+194 from CI. Also, the 
position determinations are more affected in 6 due to the propagation medium's distorting 
effects depending strongly on the elevation angle of the source. Therefore, any residual 
uncorrected ionospheric or tropospheric effects should be most visible in comparisons of the 
estimated declinations of CI relative to B2250+194 with those relative to B2252+172, since 
the angular separation of B2250+194 and B2252+172 is predominantly north-south. To 
evaluate the adequacy of the models, we plot the relative coordinates, CI — B2250+194, as 
a function of the relative coordinates, CI — B2252+172, separately in a and S. For instance, 
a straight-line slope significantly larger than unity would indicate uncorrected ionospheric 
or tropospheric effects. We plot the relative coordinates for the JPL model and for PIM in 
Figure [51 The data corrected with the JPL model were taken from Tables [2] and [3], whereas 
the data corrected with PIM are not listed. 

Least-squares fit^ to the data give slopes that are all larger than unity, reflecting the 
larger separation of B2250+194 from CI than of B2252+172 from CI. However, the slopes 
are within 1.3a of unity, except for the data in 6 for the PIM model. Here the slope is 
1.78 ± 0.20. It is smaller than a slope of ~2.5 which, if there were no ionospheric and 
tropospheric correction, would likely be expected on the basis of the ~2.5:1 ratio of the 
separations of the two sources from CI. However, this slope is significantly (~ Aa) larger 
than unity and also larger than the corresponding slope for the data corrected with the JPL 
model. The JPL model thus seems to provide a better correction for the ionosphere. 

An approximate estimate of uncorrected errors due to the ionosphere and troposphere 
for position estimates using PIM (as in Paper V for IM Peg) and an upper limit of such 
errors for position estimates using the JPL model (this paper) can be obtained from a close 
inspection of the declination data corrected with PIM (lower right panel of Figure [5]). While 
the A6 values of CI — 2252 vary over a range of 0.19 mas, the corresponding values of CI 
— 2250 vary over a range of 0.34 mas, 1.78 times larger, as also indicated by the slope. 
Therefore, the difference of the variations of 0.15 mas can likely be attributed to effects 
of the ionosphere left uncorrected by PIM. Consequently, we think that for data corrected 
with PIM the peak-to-peak contribution to the error of the position difference, CI — 2250, 
in 6 is approximately ±0.08 mas. To be conservative, we take this as a standard error. 
The corresponding contribution to the error of the position difference, CI — 2252, in 6 
approximately scales with the separation and is therefore likely to be about 0.03 mas. The 



^For simplicity we fit the data via least-squares only in the vertical direction. Fitting the data in both 
directions would not significantly alter our results. Also, since the errors in the vertical direction are similar 
to each other and the errors in the horizontal direction are similar to each other too, we did not weight the 
data in the fits and also did not plot the errors in the figures. 
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Fig. 5. — The coordinate differences (except for offsets) in a (upper row) and 8 (lower row) 
of B2250+194 - CI as a function of the coordinate differences of B2252+172 - CI for data 
for which the ionospheric effects were corrected with the JPL model (left column) and with 
PIM (right column). The data corrected with the JPL model were taken from Tables [2] and 
m The data corrected with PIM are not given in the tables. The straight lines indicate the 
least-squares fits. The slope and its statistical standard error are given for the fit in each 
figure. 
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contributions to the errors in ot that we estimated from scahng are 0.04 mas and again 0.03 
mas for the two position differences, respectively. For data corrected with the JPL model, 
we considered all these estimates for PIM as still more conservative upper limits for the 
standard deviations of these errors. 



5.3.5. Total errors 

Adding in quadrature the errors from noise, deconvolution, structure, and propagation, 
we obtain estimated standard errors in a of the position differences from B2250+194 for the 
components CI, C2, Dl, D2, and Jl of 0.06 mas and for the component Jext of 0.07 mas. The 
corresponding standard errors in 5 are 0.09 mas and 0.10 mas. The corresponding standard 
errors of the position differences from B2252+172, in both a and 5, for the components CI, 
C2, Dl, D2, and Jl are 0.05 mas and for the component Jext 0.07 mas. 



5.3.6. Comparison with empirically determined errors 

For the position differences of these components from B2250+194, our estimated stan- 
dard errors are somewhat smaller than the empirically determined standard errors which 
were mostly between the values of 0.08 and 0.10 mas in both a and 5 for these components 
(see Table |2]). The reason is the large variation of some of the data points near epoch 2001 
which increased the empirically determined errors. For the corresponding position differ- 
ences from B2252+172, all for epochs after 2001, these estimated standard errors agree well 
with the empirically determined standard errors of 0.04 to 0.06 mas in a and 5, respectively, 
for the components CI, C2, Dl, D2, and Jl and of 0.07 mas in a and 5 for the component 
Jext (see Table [3]). 



5.4. Analysis of Geodetic Group-delay Observations of 3C 454.3 and 

B2250+194 



The data from two of our sources, 3C 454.3 and B2250+194, observed in many of the 
geodetic VLBI sessions, each extending over about one day, were analyzed by one of us (L. 
P.) with the VTD/post-Solve software package. We made a weighted least-squares solution 
using all available geodetic VLBI observations of 3955 sources, including B2250 +194, 3C 
454.3, and the 212 "defining" sources from the ICRF catalogue (IFey et al.ll2009l ). made at 
157 stations from 1979 to 2008 (dubbed solution gpb_2008a). All in all we used a total of 
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6.5 million determinations of group delay from observations made simultan eously at 8.4 an d 



2.3 GHz. This solution forms our CRF. It is consistent with the ICRF2 flFev et al.l 120091 ) 



which, however, does not provide information that we need as described below. 

In particular, we estimated the coordinates of B2250+194 and 3C 454.3 for each observ- 
ing session while forcing the coordinates of each other source to be constant. We also solved 
for the positions and velocities of all stations, for polar motion and UTl parameters, and 
their rates of change, for nutation daily offsets at the middle epoch of each session and for 
numerous other parameters such as those that model clock offsets, atmosphere path delays 
in the zenith direction, and tilts of the assumed atmosphere axis of symmetry. We imposed 
the constraint that the net rotation of our estimates of source positions of the 212 defining 
sources with respect to the positions of these sources in the ICRF catalogue be zero. Such 
a choice for the constraint provides the continuity of our solution to ot her VLBI solutions. 



including the ICRF2 solution. For more details, see iPetrov et al.l (120091 ). 



The individual source positions from this solution were estimated not with respect to a 
particular reference source, but with respect to the entire ensemble of observed sources. This 
approach was possible because the geodetic VLBI sessions were designed in such a way that 
(1) a set of ~100 core sources and ~20 core antenna sites ("stations") were common to all 
sessions, and (2) the resulting estimates of source positions, station positions, and Earth ori- 
entation parameters would not be strongly correlated. The presence of sources common to all 
sessions tended to ensure the consistency of determination of the interferometer orientation 
with respect to these core sources. 

Although our list of sources has 3955 objects, the relative weight of the ~100 pre- 
dominant sources ensured that they dominate the ensemble. Therefore, we interpret the 
estimated positions of B2250+194 and 3C 454.3 as positions almost entirely with respect to 
the ensemble of the ~100 predominant extragalactic objects. 

For the solution yielding our CRF, each source was assumed to be a point source. 
However, almost all of these sources exhibit some extended structure at the milliarcsecond 
scale that may also vary over time. Since these sources were not routinely imaged and 
hence the group delays not corrected for structure effects, a position determination from 
geodetic VLBI cannot be identified with respect to a specific fiducial point in the brightness 
distribution of a source. That fiducial point is therefore in principle unidentified for any 
source with structure. In general, however, the more compact the source, the closer the 
estimated position is to the peak in the brightness distribution of the source. The source 
B2250+194 is sufficiently compact that for our purposes the position determined is effectively 
that of the brightness peak. For 3C 454.3 the fiducial point for each session is less well known 
because of the complexity of the source structure and its changes with time. Moreover, the 
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fiducial point could change due to changes of the u-v coverage. For these reasons, we do 
not rely on the 3C 454.3 position from geodetic VLBI observations, but rather use our 
dedicated VLBI observations and phase-referenced data to identify individual components 
in the brightness distribution of 3C 454.3 and to determine their positions relative to our 
two more compact extragalactic sources, B2250+194 and B2252+172, which in turn we tie 
to the CRF. 

6. Astrometric Results (2): Fit for the Position at Epoch and Proper Motion 

6.1. Fit: Components of 3C 454.3 Relative to B2250+194 and B2252+172 

We determined the position at epoch and proper motion of each of CI, C2, Dl, D2, Jl, 
and Jext relative to the brightness peaks of B2250+194 and B2252+172 with weighted least- 
squares fits for a and separately for 6. We list our results together with the weighted linear 
correlation coefficients and the weighted post-fit rms values (wrms) in Table H] (solutions #1 
and #2). We used the data corrected with the JPL ionosphere model from Tables |2] and 
El The data and the corresponding lines from the fits (solutions #1 and ^^2) are plotted in 
Figures El El El El CHI and [TTl Not surprisingly, the smallest wrms values, in both a and 
5, were obtained for the components of 3C 454.3 relative to the close reference B2252+172 
(solution #2). Our combined proper-motion estimates are given as solution 7^3. They 
were obtained as a weighted average of proper-motion estimates relative to B2250-I-194 and 
B2252+172 for the same (short) time range. 

6.2. Fit: B2252+172 Relative to B2250+194 

The position at epoch and proper motion of B2252+172 relative to B2250+194 were 
determined by first differencing the position determinations of CI relative to B2250+194 and 
CI relative to B2252+172 from Tables E] and El in the sense (CI - 2250) - (CI - 2252)1. 
We then used weighted least-squares to fit a straight line to these differences. We list the 
results also in Table HI (solution #4) and plot the data with the fit line in Figure [121 The 
relative proper motion of the two sources is zero within a small portion of la. The la upper 
limits are 11 and 24 yuas yr"'^ in a and 6, respectively. 



^Differencing the phase delays for each scan at each epoch would have given us the position of B2252+172 
relative to B2250+194 directly for each epoch and likely with a somewhat smaller uncertainty. However, our 
procedure proved to also give sufficiently accurate results for our purposes. 



~i 1 r 



"T 1 r 



m 



1.5 



1 



S 0.5 


1.5 



m 



1 



^ 0.5 




' ' ' r 
•CI - B2250+194 
o CI - B2252+172 



"T 1 T 



30±8 Aias yr-i 
20±17 /xas yr-i 



H h 



H ^ h 



H ^ \ 



H ^ h 



-27±8 /2as yr 
4± 14 /xas yr^i 

^ ^ ^ [ 



1 



^^^^ O iPffiffi^ USD 



J I L 



_l I L 



_l I L 



1998 



2000 



2002 
Year 



2004 



2006 



Fig. 6.— The coordinates of CI in 3C 454.3 relative to those of B2250+194 and B2252+172 
(except for offsets) as a function of time. The values of Aa and A5 are obtained from the 
entries in Tables [2] and [31 For discussion of apparent partial correlations between the two 
position-difference data sets, see text. 
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Fig. 7. — As in Figure El but now for C2 in 3C 454.3. 
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Fig. 8. — As in Figure El but now for Dl in 3C 454.3. 
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Fig. 9. — As in Figure El but now for D2 in 3C 454.3. 
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Fig. 11.— 



As in Figure El but now for Jext in 3C 454.3. 
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Fig. 12.— The coordinates of B2252+172 relative to those of B2250+194 (except for an 
offset in each coordinate) as determined from the differences of the position estimates, (2252- 
Cl)-(2250-Cl) given in Tables [2] and 13 
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6.3. Fit: 3C 454.3, B2250+194, and B2252+172 in the CRF 

6.3.1. 3C 454.3 

In Figure [12] we plot the position determinations of 3C 454.3 from geodetic VLBI ob- 
servations spanning almost 30 years. The statistical standard errors of the coordinates vary 
widely in part because of the different lengths of time this source was observed in the various 
sessions, and in part because the sensitivity of the VLBI systems used for the observations 
generally improved over time. 

We solved for the position and proper motion of 3C 454.3 by using weighted least-squares 
fits for a and S separately, for the whole time range of observations from 1980 to 2008. Since 
these weighted least-squares fits to the data gave xl/^ larger than unity, we again added 
a constant in quadrature to the statistical standard errors, separately for a and 6, so as to 
obtain x^/v = 1 for each coordinate. The standard errors for the individual data points so 
determined, together with the fit lines, are plotted in Figure [13] and the results are listed in 
Table [5] (solution #1). In our CRF, the proper-motion component in 5 is zero within 0.8(T. 
However, in a it is 18 ± 5 yuas yr^^, non-zero at a 3. 60" significance level. However, since 
this determination does not refer consistently to any particular component in the brightness 
distribution of the source, its significance is not clear. 



6.3.2. B 2250+ 194 

The second of our sources observed with geodetic VLBI is B2250+194. About 12 years 
of such observations yielded determinations of position and proper motion of B2250+194 
in the CRF. Via this tie and our phase-delay observations, we determined the positions of 
the components of 3C 454.3, of B2252+172, and of IM Peg in the CRF. For the source 
B2250+194, the standard errors of the position estimates also vary widely from session to 
session because of the different spans of time over which observations of this source were 
spread in the various sessions. As for 3C 454.3, we plot the coordinates for B2250-I-194 from 
the different observing sessions and their standard errors, computed as described above, 
together with the linear fit lines (solution 7^2 in Table [5]) for the whole observing period, 
from 1997 to 2008, in Figure [H] The proper motion is zero within 0.4(T. 
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Fig. 13. — The coordinates of 3C 454.3 as determined from routine geodetic VLBI obser- 
vations of up to ~4000 extragalactic sources scattered over the sky. Here and hereafter for 
plots of coordinates versus time, the straight hues give the weighted least-squares fit to the 
data points. The slopes of the lines and their statistical standard errors are indicated. These 
errors are standard errors adjusted with a constant added in quadrature to the statistical 
standard errors for each coordinate so that X^/^ = 1 (see text) . 
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The change in coordinates of B2250+194. Otherwise similar to Figure [T3l 
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6.3.3. B2252+172 

Using the position at epoch and proper motion of B2250+194 in the CRF, given in 
Table |5] (solution #2), we also obtain the equivalent values for B2252+172 in the CRF. In 
particular, we take the values from the above solution ^2 and add them to the values from 
our phase-reference observations from Table S] (solution #4) in the sense 2250 + (2252 — 
2250) and list them also in Table [5] (solution #3). The proper motion is zero within O.lo". 

Owing to a recent VLBA sensitivity upgrade, it became possible now to determine group 
delays of sources as weak as 10 mJy. Source B2252+172 was observed in two scans each 320 s 
long, with a 9-station VLBA network at epoch 2011 August 14 at 8.4 G Hz during a gap in th e 



schedule of target sources in a VLBI survey of bright infrared galaxies (ICondon et al.ll201ll ). 
The source appeared slightly resolved with a correlated flux density in a range of 11-13 mJy 
at a range of baseline projection lengths 3-200 megawavelengths. Its group-delay coordinates 
from that experiment are 22'^54'^59!597 449 ± 0.000 015, +17°33'24'.'690 49 ± 0.000 42. These 
coordinates agree with those from Table H] (solution #3) within the formers sixfold larger 
uncertainties. The corresponding proper motion for the 6.6 years between our reference 
epoch (see Table H]) and the observing date in 2011 is 13 ± 33 //as and —34 ± 64 //as, which 
gives la upper limits about threefold larger than those in Table HI 



6.4. Fit: Components of 3C 454.3 in the CRF 

We determined the position at epoch and proper motion of each of CI, C2, Dl, D2, Jl, 
and Jext in the CRF by combining the values relative to B2250-I-194 and B2252-I-172 with 
the values of the latter sources in the CRF. In particular for the positions, we take the values 
from Table |4] (solution 7^1) and add them to the values from Table [5] (solution #2) in the 
sense (3C 454.3 components — 2250) + 2250. For the proper motion values we include the 
data involving B2252+172 since they are independent of those involving B2250+194. We 
thus take all data into account, but do so only for the restricted range of epochs for which 
we have B2252+172 VLBI observations. In particular, we took the average solution #3 from 
Table m and added to it the values of B2250-I-194 in the CRF by using the solution #2 from 
Table in the sense (3C 454.3 components - 2250/2252) + 2250. We hst the resuhing 
position-at-epoch and proper-motion values for each component in Table |5] (solution #4). 
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6.4-i- Position of the Core Component, CI, of 3C 454-3 

We now discuss position determinations obtained via our two different astrometric tech- 
niques: the position of 3C 454.3 in the CRF obtained from geodetic group-delay observations 
(solution #1, Table [5]) and the position of CI in the CRF obtained through a combination 
of geodetic group-delay and our phase-delay observations (solution #4, Tabled]). We em- 
phasize that these two estimates are not expected to coincide, even in principle, because the 
former relates to some (ill defined) average over source structure and the latter to a far better 
defined component (CI) within that structure. The difference between the two position de- 
terminations (pure group delay minus combination) is 9.1 ±6.8 /zs in a and —338 ±105 yuas in 
6. While the difference in a is only 1.3cr, the difference in S is 3.2a, large enough to perhaps 
be significant. 

In this context we compare our pure group-delay p osition estimate for 3C 454.3 with 



other such estimates. For instance, the ICRF2 catalog fiFey et al.l |2009| ) provides position 
estimates of 3C 454.3 and B2250±194@ which are close to ours. It is based on essentially 
the same set of group delays, the same data editing, and the same software as we used but 
with a slightly different reduction and estimation model. The differences (ours — ICRF2) 
are: —16 ± 62 /is and —0.2 ± 3.5 fis in a and —4 ± 970 /ias and 39 ± 62 fias in S, for 
3C 454.3 and B2250±194, respectively. Here we take our estimate of an apparent proper 
motion in a of 18 ± 5 /ias yr^^ into account and propagate our estimate back 11.13 years to 
their mean epoch of 1993.95. Our estimate of the proper motion in 6 and the corresponding 
estimates for B2250±194 were small enough so that it was not necessary to consider them. 
The uncertainties are those from the ICRF2 listings only. 

Further, the latest estimate from the USNO celestial reference frame solution, crf2OO9t0, 
which uses essentially the same data set as we used only extended by another year, is different 
from ours in the sense (ours — USNO) by —1.7 ±6.1 /is and —1.4 ±1.9 /is in a and —21 ±17 
/^as and 52 ± 31 /ias in 6, for 3C 454.3 and B2250+194, respectively. Here again we take 
our estimate of an apparent proper motion in a for 3C 454.3 into account and propagate 
our estimate back 17.79 years to their mean epoch of 1987.29. The uncertainty is our proper 
motion error added in quadrature with the USNO position error. For the other differences, 
the proper motion estimates did not need to be taken into account. The uncertainties are 
those from the USNO position estimates only. 

Another estimate of the position of 3C 454.3, but not of B2250±194, was made recently 



^This source is incorrectly referred to as B2250+190 in the ICRF2 catalog. 
^http://rorf. usno.navy.mil/vlbi/ 
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with group-delay observations at 24 GHz (ILanyi et al.ll2010l ). The corresponding differences 



are —15.6 ± 8 /is and —177 ± 176 /ias. With our, statistically independent, errors added in 
quadrature, the difference in a reduces to 1.6cr and in b to 0.9o". We do not consider these 
differences to be significant. Nevertheless, note that the difference in 5 is in the direction to 
reduce the difference with our determination for CI. 

Since we compare here results from different catalogs, we point out that they are almost 
identical in their overall orientations. The differences in these orientations correspond to a 
level of only several tens of microarcseconds, which is negligible for our purposes. 

To summarize: First there is a —338 yuas (3.2o") difference in 5 between our pure group- 
delay position of 3C 454.3 and the combined group-delay phase-delay position of CI. Second, 
our pure group-delay position determinations for 3C 454.3 and, for comparison, also for 
B2250-I-194 agree within < 1.7cr with the USNO and ICRF2 position determinations. Third, 
the position of 3C 454.3 at 24 GHz in 5, while its error is large, cuts the above 3C 454.3/Cl 
discrepancy in half. We discuss these results below in § [HI 



6.J^..'2,. Limit on the Proper Motion of the Core Component CI of 3C 454-3 

Table H] (solution #1) shows that from 1998 to 2005 the proper motion of CI rela- 
tive to B2250+194 is in the southeast direction at a significance level of 3. 80" and 3Aa in 
a and S, respectively. In contrast, from 2002 to 2005 the motion of CI relative to either 
B2250-I-194 (solution not listed) or B2252-I-172 (solution 7^2) is smaller, and with an un- 
certainty corresponding to < 1.2 a, not significant, although the errors are larger. In any 
case, it appears that either B2250+194 moved to the northwest or CI indeed moved to the 
southeast, particularly during the time from 1998 to 2002. 

The cause of this apparent motion is not clear. Is it possible that the brightness peak 
of B2250+194 belongs to a jet that moved northwestward during the early time period? 
Figure E] shows that B2250-I-194 is elongated to the northwest and also slightly to the south, 
with the brightness peak located near the center of its curved structure. If there were a 
supermassive black hole located at the southern end of the structure, then the brightness 
peak would likely belong to a jet component moving away from the black hole, in this case to 
the northwest. The geodetic observations which determined the proper motion of B2250-I-194 
in the CRF show an insignificant northward motion of 21 ± 46 /zas yr~^ for the time from 
1998 to 2005 (solution not listed). The observations over the longer period from 1996 to 2008 
are collectively more sensitive, but do not indicate any motion to the northwest. Moreover, 
the geodetic observations of 3C 454.3 itself are too insensitive for a useful proper-motion 



- 34 - 



determination for the period from 1998 to 2005. 

A useful way to test this jet hypothesis would be to det ermine the position of the 



brightness peak of B2250+194 at other frequencies (see also, iKovalev et al.l |2008[ ) . The 
location of the peak at the highest radio frequencies is expected to be close to the source's 
core and the putative black hole there. For the similarly compact and elongated source, 
M81*, in the center of the nearby galaxy M81, this method did lead to the confirmation 
of the approximate location of the core of M81*. This location was earlier determined 
as the most stationary in the varying brightness distributi on of M81* relative to another 



source (the shell cente r of SN 1993 J) in the same galaxy ( iBietenholz et al.ll200ll : see also 



Bietenholz et al.ll2004j ). However, our observations at 5 and 15 GHz were not planned for 
high-precision astrometry and could not be used for this purpose. Our 5 and 15 GHz images 
appear to be, respectively, just larger and smaller versions of our 8.4 GHz images with the 
brightness peak remaining in the center of the image (Paper II), without giving any hint as 
to whether the core may be located in the south of the structure. 

As reported in Paper II, B2250+194 undergoes slight structure changes. When the 
source is modeled with an elliptical Gaussian, the major axis of the Gaussian varies between 
~0.6 and ~0.8 mas along a position angle of —11° over the time from 1997 to 2005. These 
structural changes could at least partly account for the nominal proper motion of CI relative 
to B2250+194. 

What is the la upper limit on the proper motion of CI? From the point of view of our 
IM Peg VLBI observations from 1997 to 2005 for which we use CI as a reference, the degree 
of stationarity for that period is relevant. For our longest period for the data corrected with 
JPL, from 1998 to 2005, we obtain for CI a la proper motion limit relative to B2250+194 
of 38 and 35 /las yr~^ in a and S, respectively (solution #1, Table S]). The equivalent upper 
limit of the proper motion in the CRF over that period is 46 and 56 /las yr~^ in a and S, 
respectively (solution #4, Table [5]). For GP-B (see §[1]), this is our fundamental result on the 
level of stationarity of the chosen single reference point, CI, relative to the distant universe. 

For the shorter period, from 2002 to 2005, no significant proper motion for CI was found 
above the l.ba level, neither relative to B2250 +194, nor to B2252+172, nor within the CRF. 
For this period we obtain smaller limits that are interesting from an astrophysical point of 
view. CI is stationary relative to the combination of the two reference sources within a la 
upper limit of 33 and 19 /ias yr~^ (solution 7^3 in Table [5]) and in the CRF within an upper 
limit of 39 and 30 /ias yr~^, in a and 5, respectively (solution #5, Tabled]). 
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6.4-3. Motion of the Jet Components of 3C 454-3 

What are the proper motions of the components other than CI relative to the two 
reference sources and in the CRF? Are they significantly different from that of CI and, 
moreover, can significant motion be detected of any of them relative to the distant universe? 
In Paper II we showed that C2, Dl, D2, Jl, and Jext are all moving on average away from 
CI. Their relative speeds are not necessarily constant over the more than 7 years of our 
observations. Most dramatically, Jl moves at about twice its average speed until 2002 and 
then slows down to almost zero velocity thereafter (Paper II). 

In Table H] we list the proper motions of these components relative to each of the two 
reference sources and to both combined, and in Table [5] to the CRF. For C2 we find motion 
with a larger significance than that for CI only relative to B2252+172 (solution #2) and to 
both reference sources combined {2.0a in 6, solution #3). However, we do not regard the 
significance large enough to consider the motion real. Relative to B2250+194 and to the CRF, 
C2 is stationary at about the same significance level as CI. For Jl a clear average northward 
motion of 49±9 /xas yr~^ (5.4(t) relative to B2250+194 is observed between 1998 and 2005 
(solution 7^1). The significance, however, decreases to only 1.9a, due to larger uncertainties, 
when the motion is measured in the CRF (solution #4), comparable to the significance of 
the motions in the CRF of each, CI and C2. And the motion itself decreases to (almost) zero 
within the errors for the period from 2002 to 2005 relative to either B2250+194 (solution 
not listed) or B2252+172 (solution #2) or both (solution #3) or to the CRF (solution #5) 
consistent with Jl's motion relative to CI. The largest and/or most significant motions, 
relative to either of the two reference sources or both combined or in the CRF were found 
for Dl, D2, and Jext, with, for instance, D2 having a speed in the CRF of —99 ± 25 and 
185 ± 59 /ias yr~^ in a and 6, respectively, for the short period from 2004 to 2005 (solution 
#5). 

To summarize: For the long period from 1998 to 2005, CI and C2 are stationary within 
the same small bounds relative to B2250+194 and are joined by Jl when the stationarity 
is measured (less accurately) in the CRF. The other components are stationary only within 
larger bounds or move significantly. For the short period from 2002 to 2005, CI and Jl 
are stationary within the same small bounds relative to B2252+172 and to both reference 
sources combined, and are joined by C2 when the stationarity is measured in the CRF. 
Again, the other components are stationary only within larger bounds or show significant 
proper motion. The largest significant proper motion was found for D2. 
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7. Astrometric Results (3): Analysis of Other Motions 

7.1. Fit: Parallax, proper acceleration, and orbital motion of the core 

component, CI, of 3C 454.3 

Since our computation of the motion of the guide star IM Peg includes solutions for 
parallax and proper acceleration, to help us put limits on certain systematic errors we also 
solve for these parameters for CI. In Table [HI we list the proper motion and the parallax, 
TT, obtained for CI relative to B2250+194 and separately to B2252+172 by assuming that 
both of the reference sources are infinitely distant from Earth. In addition, we solve for the 
acceleration components, jjia and /i^, for CI and also list these results in Table [61 

Our two most accurate solutions for parallax (relative to B2252+172), unsurprisingly, 
are zero (to within O.lcr), with the upper limit being ~20 /xas, corresponding to a distance, D, 
of > 50 kpc. Although other data sets might well yield more accurate parallax measurements, 
or bounds thereon, ours is one of the most, if not the most, accurate so far obtained, given the 
assumption that the reference sources are sufficiently distant to have negligible parallax. The 
proper acceleration of CI relative to B2250+194 is within lAa of zero, and is not significant. 

We also extended the number of free parameters still further and included a fit to orbital 
parameters corresponding to the period of the IM Peg binary system of 24.6 d, since such 
a fit is used in our analysis of IM Peg (Paper V). We found no indication of such orbital 
motion of CI, with each of the orbital parameters being zero within la. 

7.2. Non-linear motion of the core component, CI, of 3C 454.3 on the sky? 

Is there significant motion of the core component, CI, of 3C 454.3 on the sky that 
departs from the (linear) proper motion inherent in our fit models? Inspecting Figure [6l we 
see that the motion of CI relative to B2250+194 is somewhat correlated with Cl's motion 
relative to B2252+172 in both a and 5. To investigate possible non- linear, or in general any 
unmodeled, motion in more detail we first plot the coordinates of B2252+172 relative to 
those of B2250+194 in Figure [T5l We are using only data corrected for ionospheric effects 
with the JPL model since this model seems superior to the PIM model and since we are only 
using data taken during the period of our B2252+172 observations for which corrections 
with the JPL model were available, i.e. from 2002 to 2005. The plot shows quasi-random 
motion with wrms values of 23 /ias in a and 51 /ias in 5. In Figure [16] (upper left panel), 
we plot the positions of CI relative to those of B2250-I-194 and B2252+172. It is apparent 
that the positions of CI are covering a larger area than the positions of B2252+172 in the 
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Fig. 15. — The coordinate determinations, except for an offset, of B2252+172 relative to 
those of B2250+194 obtained by differencing the values in Tables [2] and [3] for each epoch in 
the sense (CI— 2250) — (CI— 2252). Errors are left unplotted for clarity but can be computed 
from the errors given in the above tables. The data point for the first epoch, 2002 Nov. 20, is 
indicated by a circle around the black dot. The cross in the lower left indicates the weighted 
rms of the scatter of the data points in the figure, of 23 /ias in a and 51 /ias in 5. 
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Fig. 16. — The coordinate determinations of CI, Dl, and Jext, except for an offset for eacli, 
relative to both B2250+194 and B2252+172 for each epoch (upper row). The values were 
taken from Tables [2] and [3l The lower row gives the coordinate determinations as averages 
in the sense, (Cl-2250)/2 + (Cl-2252)/2. The data points for the first epoch, 2002 Nov. 
20, are indicated in the upper panels by a circle around (i) the black dot for the motions 
relative to B2250+194 and (ii) the open circle for the motions relative to B2252+172; in 
each of the lower panels the corresponding circle is again around a black dot. The crosses 
in the lower left of each upper panel are taken from Fig. [15] and approximately indicate the 
standard errors of the data points. For the lower panels, the standard errors are assumed to 
be smaller by\/2. 
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previous figure, mainly because of a larger scatter along ot. The peak-to-peak variations are 
2.5 times larger in a but only 1.2 times larger in 5 than in the preceding figures. Also, the 
positions of CI relative to those of B2250-I-194 resemble the positions of CI relative to those 
of B2252+172. These two indications may be evidence for CI apparently moving on the sky 
along the east-west axis above the noise level which we adopt to be the wrms values from the 
B2252-I-172 versus B2250-I-194 plot and indicate by the cross in the lower left corner of the 
figure. In the lower left panel we plot the mean of each pair of position determinations (with 
-\/2 smaller cross bars) to display this apparent motion more clearly. The plotted motion is 
confined to an area not larger than ~0.2 xO.2 mas^, a small portion of the beam area but 
larger than the area of "jittery" motion in Figure [T5l Any east-west jittery motion of CI is 
consistent with our simulations (Paper II) of how CI moves east-west relative to the larger 
core region structure and relative to the 43-GHz core located 0.18±0.06 mas east of CI. This 
east-west motion simulated for CI is confined to within 0.12 mas. Our finding of possible 
motion within 0.2 mas is less precise than the result from Paper II, but in principle more 
accurate since it is measured relative to physically unrelated sources nearby on the sky. 

For comparison we show the equivalent motions of the components Dl and Jext in the 
middle and right columns of Figure [T6l As expected from our solution for the proper motion, 
Dl's motion has a linear component to the west-southwest and Jext's toward the northwest. 



8. Discussion 

8.1. Considerations for geodetic VLBI with group delays 

By comparing the positions of CI and 3C 454.3, both in the CRF, we found a possibly 
significant difference of —338 ±105 yuas (3.2o") in b while the difference in a is not significant. 
To repeat, the position of CI in the CRF was determined from the position of B2250+194 
in the CRF by the addition of the position of CI relative to that of B2250+194. Both the 
positions of 3C 454.3 and B2250+194 in the CRF were determined from geodetic VLBI 
group delays while the relative position (CI— 2250) was determined from our differential 
VLBI phase delays. In contrast to the latter, the group-delay observations were not corrected 
for structure effects. While such effects can be assumed to be insignificant for B2250+194, 
they may well be significant for 3C 454.3. Therefore, in contrast to that for B2250+194, the 
position solution for 3C 454.3 does not refer to a particular reference point in the source's 
brightness distribution. The solution should therefore be influenced by the structure of the 
source and its changes with time and frequency. In particular, a shift of the position of the 
core with decreasing freque ncy away from the foot of the jet but along the jet axis (see, 
e.g.. ISokolovsky et al.ll201ll ) could influence the results from group delays and phase delays 
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differently. 



But the source is largely oriented east-west in its brightest part, namely the core region, 
and oriented toward the northwest in its low-luminosity 10- mas long jet. On first sight 
one might therefore expect a larger discrepancy in a (see, also, iPorcasI 120091 ) and a smaller 
one, if any, in b and then with a shift to the north, not as observed to the south. If the 
10- mas jet caused the discrepancy, it would b e via some peculi ar influence. The group- 
delay position determined in 2007 at 24 GHz (ILanyi et al.l l2010l ) indeed diminishes the 5 
discrepancy in the position of CI by about half although the uncertainty of the 24-GHz 
position is large. In these observations, 3C 454.3 appears more compact than it does at 
8.4 GHz. The most easterly component at 24 GHz appears to be point-like and dominates 
the image. Low-brightness extended features ap pear up to 3 mas toward the west but 
the 10-mas jet is not visible ( ICharlot et al.l 120101 ) and s hould therefore have esse ntially no 
influence on the position determination. In this context, iMarti-Vidal et al.l (120081 ) have also 
found apparently significant discrepancies between source positions determined from geodetic 
VLBI group delays without structure corrections and those determined from differential 
VLSI phase-delays with structure corrections. All said, the cause of the discrepancy is not 
understood. 



8.2. Astrophysical implications 

8.2.1. Limit on the proper motion of the core of 3C 454-3 and the proper motions of the 

jet components 

Since the component CI is as close to stationary in the CRF as is any of our other 
five components of 3C 454.3., and could be closely related to the easternmost compact flat- 
spectrum component in 43-GHz images (Paper II), this component is likely the closest in 
our images to the expected supermassive black hole, and therefore very near the center of 
mass of the quasar. The radio emission probably originates close to the foot of the jet in the 
vicinity of this putative black hole. The jet components, C2, Dl, D2, Jl, and Jext are all 
moving away from CI (Paper II). The motions of Dl, D2, and Jext are also significant in the 
CRF at our sensitivity levels. The proper-motion values of all six components of 3C 454.3 
from Table H] (solution #5) were converted to apparent velocities and are listed in Table [71 
Our la upper limit on the proper motion of CI in the CRF of 39 and 30 /xas yr~^ in a and 
S, respectively, corresponds to a limit of 1.0 and 0.8 c. The speeds of the jet components in 
the CRF can be faster and superluminal; for instance for D2 the speed is ~ 5c. 

Similar characteristics are displayed by the superluminal quasar, 3C 345. The core was 
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found to be stationary relative to the physically unrelated quasar, NRAO 512, within 0.4c in 
the east-we st direction while t he jet components moved away from the core in this direction 
at up to 9c ( Bartel et al. 1986 1^. Since in each of these two cases the core is compact and has 
an approximately flat spectrum in the GHz frequency range, we conclude that compactness 
of the component and flatness of the spectrum are indeed, as generally assumed, indicative 
of the nearby presence of the gravitational center of the quasar. 

Component Jl is special partly because the limit on its speed over the short time interval 
from 2 002 to 2005 is sublurn inal and partly because its speed was likely much larger at earlier 
times. iPauliny-TothI (jl998[ ) found from VLBI observations at 11 GHz that a component, 
dubbed "A," moved away from the core toward the northwest from a distance of ~2.8 
mas in 1984 to a distance of ~5.0 mas in 1992 with a speed averaging 25c, but varying 
from 15c between 1984 and 1985 to 30c between 1985 and 1989 and 20c between 1989 and 



1992. This component is likely to be our Jl. IJorstad et al.l ( l2005l ) found that by 2001 this 
component, dubbed "D" in their paper, has possibly decelerated considerably to 6c. The 
strong deceleration can also be seen in our data from 1998 onward such that by 2002 to 2005 
the component remained stationary within our subluminal limits (see Table [7]). 

Larger speeds of up to 530 ± 50 /ias yr~^ were reported for jet compone nts that could 



be di scerned with higher angular resolution at 43 GHz in the CI and C2 area (IJorstad et al. 



200ll ). These authors also speculated on whether perhaps the 43-GHz core could have moved 
to the northeast by 0.1 mas in a and 0.2 mas in 8 between 1997.6 and 1998.2. While this time 
range is just before the start of our VLBI observations, our proper-motion measurements 
limit any such motion for later times. 



Recently an exceptiona lly bright optical outburst was detected in 3C 454.3 (IVillata et al. 



20061 : IVercellone et al.ll2010l ) reaching a maximum in spring 2005. It was accompanied by an 
increase of radio emission at 43 GHz from the core that started in early 2005 and reached a 
maximum in September 2005. If the outburst is associated with activity of the core, perhaps 
with the ejection of a new jet component, then the position of CI may be expected to be 
affected. However, inspection of our graphs of the temporal changes in the position of CI 
does not show any indication of a possible emerging jet component near CI, a result not 
surprising given the 2005 mid- July end of our data set. 

To repeat, any transverse linear motion on the sky found for the core of 3C 454.3 has a 
speed < 1.0c. This limit is almost as low in magnitude as the radial velocity of the quasar, 
computed from its redshift. Transverse velocities comparable in magnitude to the redshift 
velocities are not expected for quasars for a cosmological model where the dominant motion 



^We assume the same cosmological parameters as we use here for 3C 454.3. 
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is the redshift velocity due to the expansion of the universe. Any transverse velocity for 
the cores of quasars should be at least two orders of magnitude smaller than c based on 
our knowledge of peculiar motions of galaxies and galaxy clusters. The effect of acceleration 
of the solar-system barycenter toward the galactic center is also expected to be relatively 
small. However, it is plausible that a non-zero mean proper motion with respect to the CRF 
of physically unrelated sources that are separated by less than, say, a radian on the sky is 
significant because, in effect, we construct the CRF using the approximation that the position 
of the currently estimated solar-system barycenter is an inertial reference. Fortunately, the 
estimated acceleration of that barycenter toward the Galactic center has a relatively small 
effect: At its maxir num value (approx i mately appl icable for our sources), the effect is only 



about 4 A tas yr ^ (ISoyers et al. 



observed ( jTitov et al.ll201ll ) 



19981 : iTitovl |2010| ) and has recently been reported to be 



Moreover, any unexp ectedly large accelerat i on of the solar-system barycenter is less 
likely, given the study by IZakamska fc Tremaind (120051 ). who find that pulsar timing data 
(from both single and binary pulsars) are inconsistent with any unmodeled accelerations 
of the solar-system barycenter greater than ~ 4 x 10~^ cm s~^, which is only about twice 
the magnitude of the galactocentric acceleration. The more recent pulsar VLBI results of 
Deller et all (120081 ) for PSR J0437-4715 likely strengthen this upper limit. 



In the future it may be possible to search for the proper motion of the cores of quasars 
with uncertainties much smaller than c. Then it could be confirmed for the first time that 
the Hubble flow dominates the motion of quasar cores and that the velocity expected from 
the solar-system-barycenter acceleration is indeed to a high degree consistent with models. 



8.2.2. Non-linear motion of the core within the boundaries of the proper-motion limit 

Our observation of possibly significant jittery motion of the core within an area of the 
sky as small as 0.2x0.2 mas^ would be only the second time such motion has been recorded 
unambiguously for a source with core-jet structure by using as a reference for the motion 
a physically unrelated source nearby on the sky. The first source for which such motion 
was detected is the core in the core-jet structure of the nearby galaxy M81. In this case, 
the center of the expa nding shell of the nearby supernova 1993 J was used as a reference 



(IBietenholz et all 120011 1. 



What caused this apparent non-linear motion? Could such motion be indicative of 
orbital motion related to a binary black-hole system? This possibility is unlikely since the 
motion is jittery and its magnitude too large to be physically reasonable. 
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Our measurements of the jittery motion of CI along the north-south direction is within 
the noise determined from the jitter in the positions of B2250-I-194 relative to B2252+172. 
Cl's jittery motion along the east-west direction is above the noise level. In fact, the peak- 
to-peak variation is 2.5 times larger than the corresponding variation in the positions of the 
two reference sources and may be significant. If so, it is likely caused by slight brightness- 
distribution changes due to activity at the foot of the jet close to the putative supermassive 
black hole. Such changes may have influenced the flt position of the Gaussian core compo- 
nent, CI (Paper II). 

This result, if confirmed, has implications for high-precision astrometric observations 
in general. It shows that any component, even a core component, clearly identified in the 
structure of a celestial object may still move on the sky. 

8.3. The relevance for GP-B 

The goal of GP-B at launch was to measure the precession of the gyroscopes relative to 
distant inertial space with a standard error of 0.5 mas yr~^ or less in each sky coordinate. 
To be a minor contributor to the error budget, the proper motion of the guide star, IM 
Peg, was to be determined with a standard error no larger than 0.14 mas yr~^ in each sky 
coordinate. Our reference source, the quasar 3C 454.3 was shown to be stationary within 
the CRF over ~ 30 years of geodetic VLBI observations to within 0.023 and 0.009 mas yr~^ 
in a and S, respectively. More to the point, our primary reference point for GP-B, CI in 
3C 454.3, is stationary with respect to B2250-I-194 to within 0.038 and 0.035 mas yr~^ and 
within the CRF to within 0.046 and 0.056 mas yr~^ in a and 6, respectively, over the seven 
years that almost entirely cover the period of our VLBI observations in support of the GP-B 
mission and is therefore a negligible contributor to the error budget of the proper motion of 
the guide star. 

9. Conclusions 

Here we summarize our observations and data analysis, and give our conclusions: 

1. We made differential VLBI observations at 35 epochs of the quasar 3C 454.3 and the 
radio galaxy B2250-I-194 along with 12 epochs of the extragalactic, unidentified source, 
B2252-I-172, at 8.4 GHz between 1997 and 2005. With these sources we provided a reference 
frame composed of extragalactic sources nearby on the sky to IM Peg and, together with 
geodetic VLBI observations made by others we provided a (global) CRF, the latter closely 
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linked to the ICRF2, for the determination of the proper motion of the GP-B guide star IM 
Peg with respect to the distant universe. 

2. We analyzed our differential VLBI observations using phase- referenced mapping and 
phase-delay fitting in combination with a Kalman filter. 

3. Our Icr upper hmit of the proper motion of B2252-M72 relative to B2250-M94 is 
11 //as yr~^ in a and 24 //as yr~^ in 5 for the time from 2002 to 2005, identifying B2252-I-172 
unequivocally as extragalactic and providing for a highly stable reference frame of sources 
nearby on the sky to 3C 454.3. 

4. Our 1(7 upper limits of the proper motions of 3C 454.3, B2250+194, and B2252+172 in the 
CRF determined with geodetic observations, and the latter two also partly with phase-delay 
observations, are < 30 /xas yr~^ in each coordinate. 

5. Our 1(7 upper limit on the proper motion of CI, our core component of 3C 454.3, relative 
to the combination of B2250+194 and B2252+172 for the time from 2002 to 2005 is < 
35 lias yr~^ in each of a and 8, indicating that CI is highly stable with respect to two 
extragalactic sources nearby on the sky. 

6. The 1(7 upper limit on the proper motion of CI in the CRF for the time from 1998 
to 2005 is 46 and 56 //as yr~^ in a and 5, respectively.. This is our fundamental result of 
the stationarity of the reference point for the guide star IM Peg in support of the GP-B 
mission. For the shorter time from 2002 to 2005 the Icr upper limit on the proper motion 
of CI in the CRF is 39 and 30 /xas yr~^ for the two coordinates, respectively, corresponding 
to subluminal motion of < 1.0c and < 0.8c, for a and 5, respectively, for an angular-size 
distance to 3C 454.3 of 1.6 Gpc, for a flat universe with ifo=70 km s~^ Mpc~^, Q,m — 0.27, 
and Qa = 0.73. 

7. The source coordinates of CI in the CRF differ from those of 3C 454.3 determined from 
geodetic group-delay data by 131 ±98 fias in a and —338 ±105 /xas in 5, the latter coordinate 
being different at the 3.2(7 level. This difference in 5 is not understood, except possibly as a 
statistical fluke. 

8. C2 and Jl, the latter for the period from 2002 to 2005, are stationary in the CRF to 
within bounds as small as those for CI. However, Cl's flat spectrum and compactness in 
contrast to the spectra and the compactness of the other components, indicate that CI is 
closest to the putative supermassive black hole and the probable gravitational center of the 
quasar. 
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9. The jet components, Dl, D2, and Jext clearly move in the CRF. Their motions correspond 
to superluminal speeds, which for D2 is 5 c. 

10. Notwithstanding our limit on the proper motion of CI, there is evidence for its having 
jittery ~ 0.2 mas east-west motion above the noise level, likely related to jet activity in the 
vicinity of the core. This evidence is consistent with the jittery motion of CI found in Paper 
II. 

11. The la upper limit on the parallax of CI relative to B2250+194 and B2252+172 is 
20 /ias, one of the most, if not the most, accurate limit so far obtained, corresponding to 
an unsurprising lower limit of 50 kpc on its distance from Earth and demonstrating the 
sensitivity of parallax measurements with VLBI. 

12. The upper limit on the proper motion of 3C 454.3 over ~ 30 years of geodetic VLBI 
observations and of CI over ~ 8.5 years of our phase-delay VLBI observations is sufficiently 
small to meet the goal of the GP-B mission and therefore to justify use of CI as the primary 
reference point for GP-B. 
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Table 1. Characteristics of the sources 



Source 


Type 


Separation 


Flux density^ 


Redshift 


Distance ^ 










(Jy) 




(Mpc) 


3C 454.3 


quasar 






7-10 


0.859 


1610 


B2250+194 


galaxy 


-0.2 


3.6 


0.35 - 0.45 


0.28 


880 


B2252+172 


unidentified 


0.4 


1.4 


0.017 






IM Peg 


RS CVn 


-0.1 


0.7 


0.005 - 0.05 


0.0 


0.0 



^The range gives the lowest and highest flux density we measured at 8.4 GHz with 
the VLA during the course of our observations, 1997 January to 2005 July. 

^The angular diameter distance for a flat universe with 
ffn=70 km s~^ Mpc~\ and normalized density parameters, Qm = 
(IKantowski et ahlbood ). 



Hubble constant, 
0.27 and f7A=0.73 



Table 2. Separation of 3C 454.3 components from B2250+194 



Epoch 




Julian date 
2450000.0 + 




^mas) 


CI 




AS 

^mas) 






Aa 
mas) 




02" 


AS 

^mas) 






Dl 

^mas) 


AS 
(mas) 


Aa 
(mas) 


02" 

AS 
(mas) 




Aa 

^mas) 


Jl 




AS 
(mas~ 




Aa 
(mas) 


Jex 


a 


AS 
(mas' 




1998 Sep 


17 


1073 


8 


-0 


11±0 


08 





25±0 


09 


-0 


58±0 


11 





22±0 


10 












-5 


57±0 


07 





72±0 


10 


-6.34±0 


10 


2 


09±0 


OS 


1999 Mar 


13 


1251 


3 


-0 


21±0 


08 





llitO 


08 


-0 


6S±0 


12 





06±0 


09 












-5 


63±0 


07 





61±0 


09 


-6.46±0 


11 


1 


99±0 


OS 


1999 May 


15 


1314 


1 


-0 


11±0 


08 





14it0 


08 


-0 


52±0 


12 





07±0 


09 












-5 


4S±0 


07 





64it0 


09 


-6.29±0 


11 


1 


9S±0 


OS 


1999 Sep 


18 


1440 


8 


-0 


27±0 


08 





23±0 


08 


-0 


6S±0 


12 





14it0 


09 












-5 


66±0 


07 





74it0 


10 


-6.54±0 


11 


2 


14±0 


OS 


1999 Dec 


09 


1522 


5 


-0 


12±0 


08 





21it0 


08 


-0 


60±0 


12 





16±0 


09 












-5 


59±0 


07 





76±0 


10 


-6.52±0 


10 


2 


24±0 


OS 


2000 May 


15 


1680 


1 


-0 


10±0 


08 





16±0 


08 


-0 


71±0 


12 





17±0 


09 












-5 


64±0 


OS 





87±0 


08 


-6.50±0 


10 


2 


05±0 


OS 


2000 Aug 


07 


1763 


9 


-0 


13±0 


09 





lOitO 


09 


-0 


7S±0 


14 





08±0 


10 












-5 


64±0 


09 





82it0 


10 


-6.59±0 


11 


2 


17±0 


09 


2000 Nov 


05 


1854 


6 


-0 


26±0 


08 





06±0 


08 


-0 


S4±0 


13 


-0 


Olio 


09 












-5 


6S±0 


OS 





77±0 


09 


-6.53±0 


10 


1 


99±0 


OS 


2000 Nov 


06 


1855 


6 


-0 


44±0 


13 


-0 


36±0 


IS 


-1 


14±0 


16 


-0 


42±0 


IS 












-5 


93±0 


12 





35±0 


19 


-7.01±0 


14 


2 


Olio 


IS 


2001 Mar 


31 


2000 


2 





06±0 


09 





19±0 


10 


-0 


37±0 


12 





13±0 


11 












-5 


42±0 


OS 





96±0 


10 


-6.48±0 


11 


2 


34±0 


09 


2001 Jun 


29 


2090 





-0 


05±0 


08 





14it0 


08 


-0 


59±0 


14 





09±0 


09 












-5 


57±0 


09 





95±0 


09 


-6.50±0 


11 


2 


22±0 


OS 


2001 Oct 


19 


2202 


7 


-0 


12±0 


11 





02±0 


15 


-0 


70±0 


14 


-0 


05±0 


16 












-5 


65±0 


11 





88±0 


16 


-6.68±0 


13 


2 


29±0 


15 


2001 Dec 


21 


2265 


5 


-0 


21±0 


08 


-0 


06±0 


OS 


-0 


S9±0 


14 


-0 


13±0 


09 












-5 


75±0 


OS 





79±0 


09 


-6.81±0 


11 


2 


2S±0 


OS 


2002 Apr 


14 


2379 


2 


-0 


09±0 


09 





12±0 


11 


-0 


64±0 


12 





01±0 


12 












-5 


5S±0 


09 


1 


06±0 


13 


-6.65±0 


11 


2 


41±0 


11 


2002 Jul 


14 


2469 


9 


-0 


09±0 


08 





10±0 


09 


-0 


60±0 


14 





00±0 


10 












-5 


59±0 


OS 





98±0 


09 


-6.65±0 


11 


2 


34±0 


OS 


2002 Nov 


20 


2599 


6 


-0 


12±0 


09 


-0 


16±0 


10 


-0 


69±0 


13 


-0 


25±0 


11 












-5 


71±0 


09 





77±0 


11 


-6.79±0 


11 


2 


22±0 


09 


2003 Jan 


26 


2666 


4 


-0 


08±0 


09 





06±0 


10 


-0 


67±0 


13 


-0 


06±0 


11 


-1 


56±0.07 


-0.20±0.07 






-5 


65±0 


09 





95±0 


11 


-6.67±0 


11 


2 


2S±0 


09 


2003 May 


18 


2778 


1 





05±0 


08 





05±0 


08 


-0 


64±0 


12 


-0 


01±0 


09 


-1 


51±0.0S 


-0.24±0.06 






-5 


64±0 


OS 


1 


03±0 


10 


-6.59±0 


10 


2 


30+0 


OS 


2003 Sep 


08 


2891 


7 





02±0 


08 





02±0 


09 


-0 


69±0 


12 


-0 


06±0 


10 


-1 


62±0.0S 


-0.30±0.07 






-5 


60±0 


11 


1 


00±0 


11 


-6.61±0 


11 


2 


2S+0 


09 


2003 Dec 


05 


2979 


5 


-0 


01±0 


08 





11±0 


09 


-0 


73±0 


11 





01±0 


10 


-1 


70±0.09 


-0.27±0.06 






-5 


59±0 


OS 


1 


03±0 


10 


-6.63±0 


10 


2 


sSao 


OS 


2004 Mar 


06 


3071 


3 





02±0 


08 





00±0 


09 


-0 


6S±0 


12 


-0 


10±0 


10 


-1 


66±0.09 


-0.37±0.05 


-3.53±0.05 


-0.50±0.13 


-5 


56±0 


OS 





96±0 


10 


-6.62±0 


10 


2 


3S0 


OS 


2004 May 


18 


3144 


1 


-0 


03±0 


09 





09±0 


09 


-0 


71±0 


12 


-0 


06±0 


10 


-1 


69±0.07 


-0.29±0.07 


-3.59±0.05 


-0.41±0.12 


-5 


57±0 


10 


1 


05±0 


12 


-6.72±0 


11 


2 


53t0 


09 


2004 Jun 


26 


3183 





-0 


13±0 


08 


-0 


04±0 


09 


-0 


S3±0 


16 


-0 


1S±0 


10 


-1 


S6±0.07 


-0.37±0.06 


-3.69±0.05 


-0.53±0.10 


-5 


61±0 


10 





91±0 


10 


-6.73±0 


11 


2 


31+0 


OS 


2004 Dec 


11 


3351 


5 





01±0 


08 





01±0 


08 


-0 


67±0 


12 


-0 


06±0 


09 


-1 


66±0.10 


-0.44±0.05 


-3.70±0.04 


-0.39±0.10 


-5 


62±0 


09 





94±0 


09 


-6.77±0 


10 


2 


45+0 


07 


2005 Jan 


15 


3386 


4 





06±0 


08 


-0 


05±0 


OS 


-0 


62±0 


12 


-0 


10±0 


09 


-1 


66±0.10 


-0.4S±0.05 


-3.64±0.05 


-0.55±0.10 


-5 


57±0 


OS 





88±0 


09 


-6.75±0 


10 


2 


42+0 


OS 


2005 May 


28 


3519 


1 





00±0 


08 





04±0 


08 


-0 


69±0 


12 





00±0 


09 


-1 


S2±0.09 


-0.40±0.05 


-3.71±0.05 


-0.23±0.10 


-5 


63±0 


OS 


1 


02±0 


10 


-6.75±0 


10 


2 


52+0 


OS 


2005 Jul 


16 


3567 


9 


-0 


04±0 


08 





09±0 


09 


-0 


76±0 


12 





06±0 


10 


-1 


92±0.09 


-0.34±0.06 


-3.74±0.05 


-0.23±0.11 


-5 


69±0 


OS 


1 


05±0 


11 


-6.86±0 


10 


2 


60+0 


OS 



^The coordinate differences of the components of 3C 454.3 from those of B2250+194 (3C 454.3 - 2250) for each epoch for which the JPL model for the correction of the ionospheric effects could be used: 
00*^ 00"" 50^3787837 + Aa and —3° 33' 4l''067505 + A6. The coordinate differences are based on our differential measurements of CI relative to B2250+194 and on the determinations of C2, Dl, D2, Jl, 
and Jext relative to CI (Paper II). For B2250 + 194 the CRF coordinates 22^ 53"^ 7^3691736 and 19° 42' 34''628786 (solution #3, Table[5) were used. The standard errors are the statistical standard errors 
with a constant added in quadrature so that X^/^ — 1- 



Table 3. Separation of 3C 454.3 components from B2252+172 



Epoch Julian CI" C2" Dl" D2" Jl" 

date 

2450000+ Aa AS Ac AS Ac AS Aa AS Aa AS 

(mas) (mas) (mas) (mas) (mas) (mas) (mas) (mas) (mas) (mas) 



2002 Nov 20 


2599.6 


-0.11±0.05 


-0.06±0.04 


-0.68±0.06 


-0.15±0.05 










-5.70±0.05 


0.87±0.04 


-6.78±0.07 


2.32±0.07 


2003 Jan 26 


2666.4 


-0.05±0.05 


0.07±0.04 


-0.64±0.07 


-0.06±0.05 


-1.53±0.05 


-0.20±0.06 






-5.62±0.06 


0.95±0.05 


-6.64±0.07 


2.28±0.07 


2003 May 18 


2778.1 


0.04±0.05 


-0.02±0.04 


-0.65±0.03 


-0.09±0.04 


-1.53±0.07 


-0.32±0.05 






-5.65±0.01 


0.96±0.05 


-6.61±0.06 


2.23±0.06 


2003 Sep 08 


2891.7 


0.02±0.05 


0.01±0.04 


-0.69±0.04 


-0.08±0.04 


-1.62±0.06 


-0.32±0.05 






-5.61±0.08 


0.99±0.05 


-6.61±0.06 


2.27±0.07 


2003 Dec 05 


2979.5 


-0.02±0.05 


0.00±0.04 


-0.74±0.03 


-0.10±0.05 


-1.71±0.07 


-0.37±0.05 






-5.60±0.04 


0.93±0.05 


-6.64±0.06 


2.27±0.07 


2004 Mar 06 


3071.3 


-0.01±0.05 


-0.02±0.04 


-0.70±0.03 


-0.12±0.05 


-1.69±0.07 


-0.39±0.05 


-3.56±0.03 


-0.52±0.11 


-5.59±0.03 


0.93±0.05 


-6.64±0.06 


2.36±0.07 


2004 May 18 


3144.1 


-0.06±0.05 


0.06±0.05 


-0.74±0.05 


-0.08±0.05 


-1.73±0.06 


-0.32±0.05 


-3.62±0.03 


-0.43±0.10 


-5.61±0.07 


1.03±0.07 


-6.76±0.06 


2.50±0.07 


2004 Jun 26 


3183.0 


-0.09±0.05 


-0.04±0.04 


-0.79±0.11 


-0.18±0.05 


-1.82±0.05 


-0.38±0.05 


-3.65±0.02 


-0.53±0.0S 


-5.57±0.06 


0.90±0.04 


-6.69±0.06 


2.30±0.06 


2004 Dec 11 


3351.5 


0.02±0.05 


-0.04±0.04 


-0.66±0.04 


-0.11±0.04 


-1.65±0.08 


-0.49±0.05 


-3.69±0.02 


-0.45±0.0S 


-5.61±0.05 


0.89±0.05 


-6.76±0.06 


2.40±0.06 


2005 Jan 15 


3386.4 


0.04±0.05 


-0.02±0.04 


-0.64±0.04 


-0.08±0.04 


-1.69±0.09 


-0.46±0.05 


-3.66±0.03 


-0.52±0.0S 


-5.59±0.04 


0.91±0.04 


-6.77±0.06 


2.45±0.06 


2005 May 28 


3519.1 


0.00±0.05 


0.03±0.04 


-0.70±0.03 


-0.01±0.04 


-l.S3±0.08 


-0.41±0.05 


-3.72±0.03 


-0.24±0.08 


-5.64±0.03 


1.01±0.04 


-6.76±0.06 


2.51±0.06 


2005 Jul 16 


3567.9 


0.00±0.05 


0.02±0.04 


-0.73±0.03 


-0.01±0.04 


-1.89±0.07 


-0.41±0.05 


-3.71±0.02 


-0.30±0.09 


-5.66±0.03 


0.97±0.06 


-6.82±0.06 


2.52±0.06 



"As in Tablo[2]but now with B2252 + 172 as a reference. The coordinate differences (3C454.3 - 2252) are -OO"" 01™ 01!8494807 + Aa and -1° 24' 3l"l293578 + AS. For B2252 + 172 the CRF 
coordinates 22'" 54™ 59!5974430 and 17° 33' 24" 690713 (solut ion #6, Tabic [st were used. The standard errors are the statistical standard errors with a constant added in quadrature so that X^/^ — 1- 



Table 4. Relative positions at epoch and proper motions of the components of 3C 454.3, B2250+194, and B2252+172 



Component 
-reference 


(h) 


(m) 


a 




(fias yr^^) 


Pa 


wrmscK 
(^as) 


(°) 


(') 


S 


(") 




M5 

(pas yr^l) 


PS 


wrms^ 
(/^as) 


Ran} 


of 


epochs^"* 


Solution # 


CI - 2250 






50.3787837 


(19) 


30 ± 8 





60 


93 


-3 


33 


41 


.067505 


(28) 


-27 ± 8 


-0.56 


85 


1998, 


.71 - 


2005.54 


1" 


C2 - 2250 






50.3787347 


(27) 


-9 ± 11 


-0. 


16 


119 


-3 


33 


41 


.067591 


(32) 


-33 ± 9 


-0.58 


97 


1998, 


.71 - 


2005.54 


1"= 


Dl - 2250 






50.3786586 


(25) 


-126 ± 32 


-0. 


79 


74 


-3 


33 


41 


.067913 


(23) 


-89 ± 23 


-0.79 


52 


2003, 


,07 - 


2005.54 


1"= 


D2 - 2250 






50.3784787 


(14) 


-115 ± 38 


-0. 


,80 


41 


-3 


33 


41 


.067960 


(42) 


196 ± 87 


0.71 


88 


2004, 


,18 - 


2005.54 


1"= 


Jl - 2250 






50.3783932 


(19) 


-4 ± 8 


-0. 


10 


80 


-3 


33 


41 


.066488 


(33) 


49 ± 9 


0.73 


96 


1998, 


,71 - 


2005.54 


1"= 


.Text - 2250 






50.3783141 


(24) 


-51 ± 10 


-0 


72 


103 


-3 


33 


41 


.065042 


(27) 


70 ± 8 


0.87 


83 


1998, 


,71 - 


2005.54 


1'= 


CI - 2252 




-01 


01.8494857 


(14) 


20 ± 17 


0. 


,36 


45 


-1 


24 


31, 


.129358 


(16) 


4 ± 14 


0.09 


37 


2002 


,89 - 


2005.54 


2<i 


C2 - 2252 




-01 


01.8495343 


(10) 


-8 ± 14 


-0. 


19 


35 


-1 


24 


31, 


.129424 


(18) 


25 ± 15 


0.46 


42 


2002, 


.89 - 


2005.54 


2<i 


Dl - 2252 




-01 


01.8496108 


(20) 


-128 ± 26 


-0. 


,86 


59 


-1 


24 


31, 


.129789 


(20) 


-81 ± 20 


-0.80 


46 


2003, 


.07 - 


2005.54 


2^ 


D2 - 2252 




-01 


01.8497412 


(08) 


-93 ± 21 


-0. 


,89 


23 


-1 


24 


31, 


.129750 


(35) 


186 ± 72 


0.76 


72 


2004, 


,18 - 


2005.54 


2'* 


Jl - 2252 




-01 


01.8498757 


(19) 


15 ± 10 


0. 


,43 


26 


-1 


24 


31 


,128410 


(19) 


12 ± 16 


0.24 


43 


2002 


,89 - 


2005.54 


2l 


Jext — 2252 




-01 


01.8499544 


(16) 


-59 ± 21 


-0. 


,67 


55 


-1 


24 


31, 


.126911 


(25) 


96 ± 22 


0.81 


59 


2002 


,89 - 


2005.54 


2-^ 


CI - 2250/2252 










20 ± 13 


















7± 12 






2002 


.89 - 


2005.54 


3" 


C2 - 2250/2252 










-10 ± 11 


















28 ± 14 






2002 


,89 - 


2005.54 


3" 


Dl - 2250/2252 










127 ± 20 


















-84 ± 15 






2002 


,89 - 


2005.54 


3" 


D2 - 2250/2252 










-98 ± 18 


















190 ± 56 






2002 


,89 - 


2005.54 


3" 


Jl - 2250/2252 










12 ± 8 


















15 ± 14 






2002 


,89 - 


2005.54 


3° 


Jext - 2250/2252 










-59 ± 16 


















104 ± 15 






2002 


,89 - 


2005.54 


3° 


2252 - 2250 


00 


01 


52.2282694 (6) 


-2 ± 9 


-0, 


,06 


22 


-2 


09 


09, 


,938073 


(20) 


5 ± 19 


0.08 


48 


2002 


,89 - 


2005.54 


4'' 


^The eoordinates 


and proper mot 


Ions of the 


eorc 


and jet components of 3C 454.3 for the 


rcfcre 


nee ep 


och - 


- the mi 


dpoint 


of the GP-B 


mission — 


2005 Feb. 


1 (2005.08) 


, using VLBI differential 



Cn 

to 



observations of 3C 454.3 relative to B2250+194 and B2252 + 172. The parameters are derived from weighted least-squares linear fits. The uncertainties are standard errors derived from the fit 
on the basis of the statistical standard errors of the individual measurements added in quadrature to constants so that X^/^-' — 1- We searched for any sign that a possible correlation betw^een 
consecutive data points (see, e.g., Fig[6) could render our uncertainties too small by repeating the weighted least squares fit for even and odd numbered data points separately. However, we did 
not find any such sign and therefore think that any possible correlation would only have a minor effect on our uncertainty estimates. Also listed are the weighted linear correlation coefficients, p, 
from the fits, and the weighted rms values, wrms, for the post-fit residuals, for the two coordinates. All data are corrected for the effects of the ionosphere by using a model implemented in the 
AIPS software and based on GPS data (JPL, see text). The earliest date at which the model could be used is 1998 Sept. 17 (1998.71). For Dl and D2 the time range starts with the first epoch 
at which a component could be identified in the brightness distribution, that is, 2003 Jan. 26 (2003.07) and 2004 Mar. 6 (2004.18), respectively. 

'^The range of epochs from which the data were taken for the solution. The range 1998.71 — 2005.54 is the time range from 1998 Sep. 17 to 2005 July 16 of our GP-B VLBI observations for 
which the ionospheric model, JPL, could be used. The epoch 2002.89 refers to 2002 Nov. 20 when we started to include B2252-|-172 in our observations. The other time ranges refer to the total 
ranges for which data were collected for the respective component and source. The reference epoch 2005.09 is 2005 Feb. 1, the midpoint of the time period through which data were taken on the 
GP-B spacecraft. 

'^Thc difi'erences from B2250-I-194 in the coordinates of the positions and in the corresponding components of the proper motion for each of 3C 454. 3's six core and jet components. The time 
range starts for all components with the first epoch of our VLBI observations, except for components Dl and D2, where it starts with the first epoch at which a component could be identified in 
the brightness distribution. The epochs of the time range arc 1998 Sept. 17 (1998.71), 2003 Jan. 26 (2003.07), and 2005 Jul. 16 (2005.54). 

*^As in c for solution #1 but now relative to B2252-K172. The time range starts with the date of the first B2252-K172 observations, 2002 Nov. 20 (2002.89). 

"^As for solution #1 but now starting not earlier than at epoch 2002.89, and including the values for B2252-|-172. In particular, for each component, we took the weighted average of (a) the 
proper motion from a solution (3C 454.3-2250) for the short time range (not listed), and (b) the proper motion of solution #2 (3C 454.3-2252). The errors from (a) and (b) were added in 
quadrature. 

^Thc coordinates of B2252-|-172 relative to those of B2250-I-194 and their changes with time. They were derived from the data (2252 — Cl) — (2250 — CI) determined for each epoch. 



Table 5. The positions at epoch and proper motions in the celestial reference frame (CRF) 



Ooiup oiicnt-rcicrcricc 








(jias yr ) 


Pa 




I ) 


I ) 


PS 
(pas yr) 


PS 


wnns ^ 


Jr\.3,l"lgC OI 


b 

cpo CIIS 


C 1 f M. 

oolULion 


3C 454.3 


22 


53 


57.7479664 (60) 


18 ± 5 


0.11 


742 16 


08 


53.560943 (80) 


4 ± 5 


0.03 


918 


1980.28 - 


2008.47 


1 


B2250+194 


22 


53 


7.3691736 (24) 


3 ± 10 


0.05 


179 19 


42 


34.628786 (62) 


-7 ± 20 


-0.06 


345 


1996.01 - 


2008.36 


2 


B2252+172 


22 


54 


59.5974430 (25) 


1 ± 13 




17 


33 


24.690713 (65) 


-2 ± 28 






2002.89 - 


2005.54 


3" 


CI (50) 


22 


53 


57.7479573 (31) 


33 ± 13 




16 


08 


53.561281 (68) 


-35 ± 21 






1998.71 - 


2005.54 


4d 


C2 (50) 


22 


53 


57.7479083 (36) 


-7± 15 




16 


08 


53.561195 (70) 


-40 ± 22 






1998.71 - 


2005.54 


4d 


Dl (50) 


22 


53 


57.7478322 (35) 


-123 ± 34 




16 


08 


53.560873 (66) 


-96 ± 30 






2003.07 - 


2005.54 


4d 


D2 (50) 


22 


53 


57.7476523 (28) 


-112 ± 39 




16 


08 


53.560826 (75) 


188 ± 89 






2004.18 - 


2005.54 


4d 


Jl (50) 


22 


53 


57.7475668 (31) 


-1 ± 13 




16 


08 


53.562298 (70) 


42 ± 22 






1998.71 - 


2005.54 


4<1 


Jext (50) 


22 


53 


57.7474877 (34) 


-49 ± 14 




16 


08 


53.563744 (68) 


63 ± 21 






1998.71 - 


2005.54 


4<i 


CI (50/52) 








22 ± 17 










4 ± 26 






2002.89 - 


2005.54 


5= 


C2 (50/52) 








-8 ± 15 










25 ± 27 






2002.89 - 


2005.54 


5= 


Dl (50/52) 








-125 ± 23 










-92 ± 27 






2003.07 - 


2005.54 


5" 


D2 (50/52) 








-99 ± 25 










185 ± 59 






2004.18 - 


2005.54 


5" 


Jl (50/52) 








13 ± 14 










10 ± 24 






2002.89 - 


2005.54 


5° 


Jext (50/52) 








-58 ± 25 










100 ± 26 






2002.89 - 


2005.54 


5° 



^Thc source coordinates, o: and S, of 3C 454.3, B2250+194, and B2252+172 at the reference epoch 2005.08, and their changes with time, /.^q and fi§ respectively, in the CRF. The parameters are 
derived from weighted least-squares linear fits. The uncertainties arc standard errors derived from the fit on the basis of the statistical standard errors of the individual measurements added in quadi^^ure 
to constants so that X^/^ — Also listed arc the weighted linear correlation coefficients, p, from the fits, as well as the weighted rms values, wrms, for the post fit residuals, for the two coordinates. 

'^The range of epochs as in Table 

"^The coordinates of B2252-|-172 and their changes with time derived by adding the solutions from #2 to those of #4 in Table g] The errors were added in quadrature. 

"^The coordinates and proper motions of the components of 3C 454.3 in the CRF. Wc added the values from solution #1 of Table [4] to the values from solution #2 of this table in the sense (3C 454.3 
comp. — 2250) -j-2250, and added the errors in quadrature. 

^As for solutions #4 but now only for the proper motion, and at epochs not earlier than 2002.89, and with the inclusion of the values for B2252 + 172. In particular, we added the values from solution 
#3 of Table to the values from solution #2 of this table in the sense {3C 454.3 comp. — 2250/2252) + 2250 and added the errors in quadrature. 



- 54 - 



Table 6. Proper motion, parallax, and proper acceleration' 



Source-reference 




Ai5 


vr 








(//as yr^i) 


{^as yr~i) 


(/ias) 


(/ias yr"2) 


(/ias yr"2) 


Cl-2250 


29± 8 


-27 ± 8 


40 ±20 






Cl-2252 


21±18 


4± 14 


-2 ±17 






Cl-2250 


33±26 


11 ±27 


38 ±20 


2 ± 9 


13± 9 


Cl-2252 


3±44 


13 ±37 


-2 ±18 


-21 ±47 


11± 40 



'^Parameter estimates for component CI of 3C 454.3 relative to B2250±194 and, sep- 
arately, to B2252±172. Uncertainties are statistical standard errors derived from the 
weighted least-squares fit, scaled to x^/v = 1. The reference epoch is 2005 February 
1. 



Table 7. Velocity of the components of 3C 454.3 in the CRF' 



Component 






Range 


of epochs^ 


CI 


0.6 ±0.4 


0.1 ±0.7 


2002.89 


- 2005.54 


C2 


-0.2 ±0.4 


0.6 ±0.7 


2002.89 


- 2005.54 


Dl 


-3.2 ±0.6 


2.3 ±0.7 


2003.07 


- 2005.54 


D2 


2.5 ±0.6 


4.7 ± 1.5 


2004.18 


- 2005.54 


Jl 


0.3 ±0.4 


0.3 ±0.6 


2002.89 


- 2005.54 


Jext 


-1.5 ±0.6 


2.5 ±0.7 


2002.89 


- 2005.54 



^Velocity in a and 5 in units of the speed of light for angular 
velocities in Table [5] (solution #5). The velocities were computed 
for an angular diameter distance of 3C 454.3 of 1.6 Gpc. 

'^Same as in Table |H 



